Reconstruction of modified gravity with ghost dark energy models 
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Abstract 

In this work, we reconstruct the f(R) modified gravity for different ghost and generalized- 
ghost dark energy (DE) models in FRW flat universe, which describe the accelerated expansion 
of the universe. The equation of state and deceleration parameter of reconstructed f(R) gravity 
have been calculated. The equation of state and deceleration parameter of reconstructed f(R)- 
ghost/generalized-ghost DE, have been calculated. We show that the corresponding f{R) gravity 
of ghost/generalized-ghost DE model can behave like phantom or quintessence. Also the transition 
between deceleration to acceleration regime is indicated by deceleration parameter diagram for 
reconstructed f(R) generalized-ghost DE model. 
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I. INTRODUCTION 



The acceleration expansion of the present universe has been accepted by many observa- 
tional evidence Although the dark energy (DE) scenario with a negative pressure can 
drive this expansion, the modified gravity, so called f(R) gravity, has been introduced as 
a self consistent and nontrivial alternative way to explain early time inflation and present 
acceleration expansion. In f(R) gravity, the Einstein-Hilbert action is modified and gener- 
alized by incorporating a new phenomenological function of Ricci scalar, f{R). 

On the other hand the dynamical DE models are proposed in order to solve DE problems 
(cosmic coincidence and fine tuning problems). Recently, the Veneziano ghost DE has 
been attracted a deal of attention in the dynamical DE category. The Veneziano ghost 

n 

is proposed to solve the U(l) problem in low-energy effective theory of QCD [2] and has 
no contribution in the flat Minkowski spacetime. In curved spacetime, however, it makes 
a small energy density proportional to Aq CD H, where Aqcd is QCD mass scale and H is 
Hubbel parameter. This small vacuum energy density can be considered as a driver engine 
for evolution of the universe. It should be mentioned that the Veneziano ghost DE model 
does not violate unitarity, causality, gauge invariance etc js-5]. In fact the description of DE 
in terms of the Veneziano ghost is just a matter of convenience to describe very complicated 
infrared dynamics of strongly coupled QCD. In other words, the veneziano ghost is not a 
new physical propagating degree of freedom. The same dynamics can be described without 
using the ghost, with some other approaches (e.g. direct lattice simulations). Also note 
that this model is totally arisen from standard model and general relativity. Therefore one 
needs not to introduce any new parameter or new degree of freedom and this fact is the 
most advantages of ghost DE. With Aqcd ~ lOOMev and H ~ 10~ 33 ef , the right order of 
observed DE density can be given by ghost DE. This numerical coincidence also shows that 



this model gets ride of fine tuning prob 
model with energy density as p = aH 



em 
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5|. Many authors have already suggested DE 



3] and p = aH + 0(H 2 ) [81]. Generally it is very 
difficult to accept such a power like behavior as QCD is a theory with a mass gap determined 
by the scale ~ 100 Mev. The power like scaling ~ H is due to the complicated topological 
structure of strongly coupled QCD, not related to the physical massive propagating degrees 
of freedom. In fact the linear in Hubble constant "if" scaling is not a baseless assumption, 
but rather has a strong theoretical support tested in a number of models, and tested in 
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the lattice QCD simulations 



Recently the ghost-scalar field models have been 



investigated 9j. Further more ghost DE model has been fitted with current observational 
data including SNIa, BAO, CMB, BBN Q,[lO|. 

On the other hands, a consistent f(R) function in modified gravity context, has been 



reconstructed by many dynamical DE models [llMl3l] . It is worthwhile to mention that the 
modified f(R) gravity model considered as an effective description of the underlying theory 
of DE, and considering the ghost DE scenario as pointing in the direction of the underlying 
theory of DE, it is interesting to study how the f(R) gravity can describe the ghost DE 
density as effective theory of DE models. This motivated us to establish the different models 
of f(R) gravity according to the ghost and generalized-ghost DE model. In this paper we 
want to reconstruct a consistent modified gravity so that it gives the cosmological evolution 
of ghost DE model. 



II. GENERAL FORMALISM OF F(R) GRAVITY 



The general action of f(R) gravity in Planck mass unit in which h = c = 8nG = 1, is 
given by [14j 



S 



-g d 4 x 



R + f(R) 



(i) 



where g, R and C m are the determinant of metric g^, Ricci scalar and matter contribution of 
the action, respectively. Here G is gravitational constant and f(R) is an arbitrary function 
of R. In metric formalism, Where the action is varied only by the metric, by taking the 
variation of action with respect to the metric g^ u , the following field equation is obtained 



Rju> 2 -^9 



-g^f(R) - R^f'(R) + (V„V„ - g^a)f(R) 



(2) 



Here the prime denotes a derivative with respect to R and R^ and T^/ are the Ricci tensor 
and the energy-momentum tensor of the matter, respectively. 

In a spatially flat-FRW background universe, with following line element 

ds 2 = -dt 2 + a 2 it) [dr 2 + r 2 d6 2 + r 2 sin 2 (#)d</> 2 ] , (3) 



by taking the matter in the prefect fluid form, in which X^"^ = diagf— p m , p m , p m , p m ), the 



field equation (j2J) reduces to following modified Friedmann equations 

1 



H 



-,Peff, 



(4) 



2H + 3H 2 = -p eff . 
Here H = a/a is the Hubble parameter and 



Peff 



Peff 



-f(R) + 3(H + H 2 )f'(R) - 18(AH 2 H + HH)f'(R) 



-f(R)-(H + 3H 2 )f(R) 



+ Pr, 



(5) 



(6) 



+6(8H 2 H + 6HH + AH 2 + H)f"(R) + 36(H + AHH) 2 f"\R) 
The Ricci scalar is obtained as 

R = 6(H + 2H 2 ), 

where the dot denotes a derivative with respect to cosmic time t. 
The energy conservation law gives 



+ Pr, 



(7) 



(8) 



Peff + 3Hp eff (l + Wtot) = 0, 



where 



Wtot 



Peff 
Peff 



(9) 



(10) 



is the total equation of state (EoS) parameter. In the absence of any matter field, L F(R) 
dominated universe', the EoS parameter has only the gravitational contribution as 



w R 



Hf\R) + 3(3HH-4H 2 H + 4H 2 + H)f%R) + 18(H + 4HH) 2 f''(R) 
f(R) —6(H + H 2 )f'(R) + 36[AH 2 H + HH)f'(R) ' 



In this case the first modified Friedmann equation (J3J) yields 

3# 2 = PR, 



(12) 



where pR is the gravitational contribution of energy density (E]) for p m = 0. 

Taking time derivative of Eq. (j!2p and using continuity equation one can find a new 
relation for EoS parameter as 



w R = -1 - 



2H 
3IP' 



(13) 



For H > 0, we have wr < — 1, indicating the phantom DE dominated universe. In the case 
of H < 0, we have wr > —1, representing the quintessence DE dominated universe. The 
deceleration parameter in the F(R) dominated universe is obtained as 

H 3w R + 1 



aa 



-1 



H 2 



(14) 



In the case of wr < —1/3, we have q < indicating the accelerated expansion of the universe 

and for wr > —1/3, we have q > 0, representing the decelerated expansion phase. 

In the theories of modified gravity there are two class of scale factor that have been usually 



used in the literature 



16] . The first one has been given by 



Using Eqs. 



a(t) = a (t s -t)~ h , t<t s , h>0. 
and (lT5"j) one can obtain 

H 



h 



U-t 



hR 



12/i + 6 



1/2 



H = H 2 /h, 



H 



2H 3 



(15) 



(16) 



Here a big rip singularity will be happened at t — t s [171 ] . It it is easy to see that for the 
first class of scale factor given by ( fl5|) we have H > 0. Hence from (Tl3|) . one can conclude 
that for this class of scale factor wr < — 1, indicating the phantom DE dominated universe. 
This scale factor is usually so-called phantom scale factor. 
The second class of scale factor bas been defined as 



a(t) = a t h , h > 0, 



The Hubble parameter, in this case, is obtained as 



t 



hR 



1/2 



H 



-H 2 /h, H 



2H 3 



(17) 



(18) 



12/i-6 

From (fl3|) one can see that H = —H 2 /h < represents the quintessence dominated 
universe. Hence, the model ffl7|) is so-called the quintessence scale factor. 



In sections (3) and (4), by using the above classes of scale factor, we reconstruct the 
different f(R)- gravities according to the ghost and modified ghost DE models. 



III. F(R) RECONSTRUCTION FROM GHOST DE 



The energy density of ghost DE is given by 



Pa = otH 



(19) 



where a; is a constant of model. 
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A. phantom scale factor 



For the first class of scale factor (fT5l) and using (Tl6|) . the energy density of ghost DE in 
( TrU]) is written as 

/ hR 

" A = a Vm + 6 < 20) 

Equating (|20|) and ([6]) by using (fl6l) gives the following solution 

/(i?) = + \_R m - - aj p R^ (21) 

where 

h + 3 ± V/i 2 - lOh + 1 , s 

m± = - (22) 

and 



2 /12/1 + 6 . , 

-*=bV — s— (23) 

Here A± are the constants of integration and can be obtained from the initial conditions. In 
order to generating accelerated expansion at the present time, one can assume that f(R) is 
a small constant at the present time as follows 

f(R ) = -2R , and f (R ) ~ 0, (24) 



where Ro ~ (10 33 ev) 2 is the current curvature [18]. Using the boundary condition (1241) . the 
constants A± in ( 12T|) are obtained as 

4v/i^m_ + o7p(l - 2m_) 

A + = ( 25 ) 

2R 2 (m + — m_) 

and 



A = 4yfipm + + a7 p (l - 2m + ) 
2i?™~ 2 (m + — m_) 

Inserting (|2T|1 into (fTTj) and using (fl6l) obtains the EoS parameter of ghost f(r) gravity as 

u, s = -l-|, (27) 

which describe the phantom phase of accelerated universe. It must be mention that, from 
Eq( l32l) . the reconstruction can be performed provided that 

< h < (5 - 2^6) or h>(5 + 2V6) (28) 
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The EoS parameter (|27|) .shows that the reconstructed f(R) gravity from the ghost DE can 
cross the phantom divide, while in the scenario of DE, the ghost model can cross the phantom 
divide if the interaction between dark matter and DE is included 19]. 



Using fTST] HI]) , the deceleration parameter is constant and calculated as 

« = "I - i (29) 

It shows that, in this case, for all values of h in the range ( I28p . our universe is expanded 
under an accelerating phase q > 0. 



B. Quintessence scale factor 

Now we reconstruct ghost f(R) gravity for the second class scale factor in (|17|) . By using 
()18p . the energy density of ghost DE in ffl~9l) can be written as follows 



PA=a Vi^b (30) 

Equating ( 1301) and (jHJ) by using ( {TBI) gives the following solution 

f{R) = X+R m+ + X-R m ~ - aj q R* (31) 

where 

3 - h ± Vh 2 + Wh + 1 , s 

m± = (32) 

and 



2 12h-6 

^ = 9y^r- (33) 

Using the boundary condition (|24l) . the constants A± in (J3T1) can be obtained as 

A = 4y/^m_ + aj q (l - 2m_) 
2 J R™ + "^(m+ -m_) 



and 



A = 4y/^m + + a7 9 (l - 2m + ) 



2Rq 2 (m + — m_ 



Inserting (13T|) into (11 ip and using (1181) results the EoS parameter of ghost f(r) gravity as 

w R =-l + ^, h>l, (36) 
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which corresponds to the quintessence regime, i.e., — 1 < wr < —1/3. In this case, the 
deceleration parameter is constant and obtained as 

9 =-l + i (37) 

In this manner we find that in quintessence regime with h > 1, the parameter q < 0, while 
for 1/2 < h < 1, we have g > 0, which reveal a deceleration expansion at early time. 
Therefore by this reconstruction, both deceleration and acceleration have been permitted, 
separately for various the parameter h. 



IV. F(R) RECONSTRUCTION FROM GENERALIZED-GHOST DE 

Here we reconstruct the f{R) gravity for generalized-ghost DE. The energy density of 
phenomenological model is given by 

p A = aH + 8H 2 (38) 



where a and 8 are the constants of the model. This model was first proposed in [20| to get 

n 

an accelerating universe. For othe motivation to consider this form see |21|. When 8 0, 
we recovered the model discussed in Sec. IIIII and for a — > 0, this model give the holographic 
DE model with Hubbel scale as IR-cutoff. Also the subleading term H 2 in the ghodt DE 
model might play a crucial role in the early evolution of the universe jlO]. In fact in this 
model we have a 'XghostDE' model where 'X' denotes the second term in ( |38l) . 



A. Phantom scale factor 

Same as previous section, first we choose the phantom scale factor ( |T5l) . Therefor the 
energy density of generalized-ghost DE, by using f|T6|) . is rewritten as 



hR n hR , . 

" A = a Vl2AT6 +/! 12ftT6 (39) 

Equating (|39|) and ([6]) by using (|T6j) results the following solution for f(R) 

f(R) = \+R m+ + X-R m - - a lp R^ - ^R (40) 
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where parameters m± and 7 P are the same as Eqs. (1221) and (|23|) . respectively. Applying 
the boundary conditions in (124]) . the constants A± in PUI) are determined as 



2^V^(1 - m_) + 3a 7p (l - 2m_) + 12^m_ 

+ ~~ ^rn i 4i ) 

6i? 2 (m+ — m_) 



A 2/3y^(l - m+) + 3a 7p (l - 2m + ) + 12Vfiom + 

6i?™~ 2 (m + — m_) 

Substituting (J1D|) into (TTTT) and using ffTB"]) yields the EoS parameter of generalized-ghost 
f(R) gravity as 

2 / R0 \ 

WR = ^ ~ 3h V + YlaE + 6aHh~i + 0R J " (43) 

The above EoS parameter is time-dependent. Here in contrast with constant EoS parameter 
in (12"?]) . The dynamical behavior of EoS parameter is given and we see that at infinity it 
merge to a constant value (l2~Tj) of Sec. IIHI 

In this case, the time varying deceleration parameter is obtained as 

! i R£ , . 

q h 12ahH + QaH + /3hR' [ ' 

Same as the previous section, this give us the accelerating expanding universe for all values 
of h in the range ( 128]) . 



B. Quintessence scale factor 

We now choose the quintessence scale factor f[T7|) . In this case the energy density of 
generalized-ghost DE, by using f JT8|) . is rewritten as follows 



pA = a y 127^6 + /? 12/^6 (45) 
Equating (145]) and ([6]) by using (fl8j) results the following solution for f(R) 

f(R) = X+R m+ + X-R m ~ - a lq R^ - ^R (46) 

where parameters m± and 7 P are the same as Eqs. (1321) and (|33|) . respectively. Applying 
the boundary conditions in (|24p . the constants A± in fj40l) are determined as 



2/3 v ^(l-m_) + 3a7 g (l-2m_) + 12 v / ^m_ 

+ ~~ t 47 J 

6i? 2 (m+ — m_) 

9 



0.5^ 








0.4: 








q 0.3: 








0.2: 








0.1: 








0; 

-0.1: 


0.2 0>- 


vQ,6 0.8 1 1.2 1.4 1.6 1.8 
a 


2 



FIG. 1: The evolution of deceleration parameter, q, versus a in generalized-ghost DE model for 
parameters: (h = 1.3, a = 2, /3 = 0.3). 



2/3v^o(l - m+) + 3«7 g (l - 2m + ) + 12v^m+ 



(48) 



m_ 



Substituting (H6|) into (TTTj) and using (II5]I obtains the EoS parameter of generalized-ghost 
/(i?) gravity as 

2 / \ 
U,fl = - 1 + 3fe( 1+ 12aff-6a^-i + /3fl j' (49) 
which can represent the quintessence regime of accelerated universe. Here we see that 
the EoS parameter is dynamical, in contrast with constant EoS parameter in (15511 . The 
dynamical behavior of EoS parameter is achieved and we see that at infinity it merge to a 
constant value (155]) of Sec. IIHI 

The deceleration parameter, in this case, is time varying and obtained as 



q = -1 + T + 



R(3 



(50) 



h YlahH + 6aH + (3hR 

Here in contrast with (157)) in previous section, the deceleration parameter evolves with time. 
In Fig. [H for h > 1, an early time transition from the deceleration to acceleration phases is 
demonstrated. 
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V. CONCLUSION 



The reconstruction of f(R) modified gravity for different ghost and generalized-ghost DE 
models in FRW flat universe has been investigated. The reconstruction was performed by 
considering two classes of scale factors, containing i) phantom scale factor, a = ao(t s — t) h 
and ii) quintessence scale factor, a = aot h . The equation of state and deceleration parameter 
of reconstructed /(-R)-ghost /generalized-ghost, have been calculated. We showed that the 
corresponding f(R) gravity of ghost /generalized ghost DE model can behave like phantom 
or quintessence. In /(i?)-generalized ghost case, the EoS parameter can vary with time. 
Also, the transition between deceleration to acceleration regime is happened at early time. 
These behaviors are in contrast with Pure /(-R)-ghost and /(i?)-holographic DE [12|, Il3 |. 
Therefore the generalized-ghost DE model give us a stronger view of the universe in f(R) 
reconstruction point of view. 
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